Super-Resolution at the Shot-Noise Limit 
with Coherent States and Photon-Number-Resolving Detectors 
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There has been much recent interest in quantum optical interferometry for applications to 
metrology, sub- wavelength imaging, and remote sensing, such as in quantum laser radar (LADAR). 
For quantum LADAR, atmospheric absorption rapidly degrades any quantum state of light, so 
that for high-photon loss the optimal strategy is to transmit coherent states of light, which suffer 
no worse loss than the Beer law for classical optical attenuation, and which provides sensitivity 
at the shot-noise limit. This approach leaves open the question — what is the optimal detection 
scheme for such states in order to provide the best possible resolution? We show that coherent 
light coupled with photon number resolving detectors can provide a super-resolution much below 
the Rayleigh diffraction limit, with sensitivity no worse than shot-noise in terms of the detected 
photon power. 
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Ever since the work of Boto et al. in 2000, it has been 
realized that by exploiting quantum states of light, such 
as NOON states, it is possible to beat the Rayleigh diffrac- 
tion limit in imaging and lithography (super-resolution) 
while also beating the shot-noise limit in phase estima- 
tion (super-sensitivity) [3, 0, However such quan- 
tum states of light are very susceptible to photon loss 
@j Hi @| • Recent work has shown that in the presence of 
high loss, the optimal phase sensitivity achieved is always 
Aip = \j\ffi, where n is the average number of photons 
to arrive at the detector 0, @]. These results suggest 
that, given the difficulty in making quantum states of 
light, as well as their susceptibility to loss, that the most 
reasonable strategy for quantum LADAR is to transmit 
coherent states of light, \a) = e~H 2 / 2 EfeLo( a Vv^0l fe > 
where |fc) is a fc-photon Fock state to mitigate a super- 
Beer's law in loss and maximize sensitivity [3]. It is well 
known that such an approach can only ever achieve at 
best shot- noise limited sensitivity However, such a 
conclusion leaves open the question as to what is the 
best detection strategy to optimize resolution. Recent 
experimental results have indicated that such a coherent- 
state strategy can still be super-resolving, provided a 
quantum detection scheme is employed [10(. In this pa- 
per we derive an optimal quantum detection scheme, 
that is super-resolving, and which can be implemented 
with photon number resolving detectors. Our proposed 
scheme exploits coherent states of light, is shot- noise lim- 
ited in sensitivity, and can resolve features by an arbi- 
trary amount below the Rayleigh diffraction limit. Our 
scheme would have applications to quantum optical re- 
mote sensing, metrology, and imaging. In Fig. 1 we illus- 
trate schematically a two-mode interferometric quantum 
LADAR scheme. The source at the left is assumed to 
contain a laser producing a coherent state \cx)a m upper 
mode A with vacuum in lower mode B, which is illus- 




FIG. 1: Here we indicate the Mach-Zehnder interferometer 
used in the calculations. The coherent state is incident in 
mode A and vacuum in mode B at the left at Line I. Af- 
ter the first beam splitter transformation we have the two- 
mode coherent state of Eq. (1), as indicated at Line II. After 
the phase shifter <p this state becomes the two-mode coherent 
state of Eq. (2). At Line III we also implement the detection 
schemes corresponding to the operators Nab , vab , and /Lab ■ 
Finally after the final beam splitter, we implement the parity 
operator Ha detection of Eq. (13) in the upper mode. 



trated at line I in Fig. 1. The state is incident on a 50-50 
beam splitter (BS), which mixes this coherent state with 
the vacuum state \0)b in lower mode B. The output of 
such a mixing is computed by the well known beam split- 
ter transformation, and is the state, 
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where the notation \p,q) = \p) Bi n = M 2 , and with- 
out loss of generality, a BS phase factor of i = e J7r/ " 2 has 
been dropped for clarity. (This physically corresponds to 
a symmetric 50-50 BS with a quarter-wave plate in one 
of the output ports.) This is the state at line II in Fig. 
1. As shown in Fig. 1 the upper mode A imparts a phase 
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shift ip on this state, yielding, 
W)ab = \ae^/V2,a/V2) 
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where, without loss of generality, we take a = y/E, which 
amounts to setting an irrelevant overall phase to zero. 
This is the state at the line III in Fig. 1. In the coher- 
ent state basis, this state is obviously separable. How- 
ever it is less obviously separable in the number basis, 
and it is easy to see that the double sum contains the 
NOON states \N ::0) v = e lNip \N,0} + \0,N) as well as 
MM' states \M::M')v = e m ^\M,M') + e iM ' v \M' , M), 
both of which are known to be N = M — M' fold super- 
resolving. Indeed this observation suggests a strategy, 
similar to that employed by Resch et at, of projecting 
the state \t/j) ab of Eq. (2) onto the NOON basis through 
the implementation of the detection operator Nab — 
\N,0)(0,N\ + \0,N)(N,0\ 10]. This results in an expecta- 
tion value of ab(iP\Nab\4>)ab = (n/2) N e~ n 2 cos Nip/ N\ 
that is clearly iV-fold super-resolving. However the fac- 
tor of (n/2) N e~ n /N\ indicates that the visibility of this 
expectation value, as a function of ip, is much less than 
unity. Let us optimize this visibility by tuning the return 
power n. The maximal visibility occurs when n = N 
or when the return power is equal to the desired super- 
resolution factor. For example, with n = N = 10 we 
achieve 10-fold super-resolution in ab("4>\Nab\iI>}ab with 
a visibility of about 0.024%. In a similar fashion we may 
now estimate the minimal phase sensitivity variance via 
the usual Gaussian error propagation formula Q , 
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where the factor of 1/N 2 would provide Heisenberg lim- 
ited sensitivity, as it would be in the case of a pure NOON 
state with this detection scheme Q , if it were not for the 
Poisson weight factors inherited from the coherent state. 
This expression, Eq.(3), can be minimized by inspection, 
again by taking N = n and simultaneously choosing ip 
(which can be done by putting a phase shifter in arm B) 
such that Nip = n/2, yielding, 
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where we have used Stirling's approximation, and where 
again the Heisenberg limit behavior of 1/n 2 appears in 
tantalizing fashion. However, the exact expression has a 
minimum at n = 2 of A</? 2 y|n=2 = 1.85, and hence for 
all values of n always does worse than the shot-noise and 
Heisenberg limits of A(,5g NL = 1/n > A^jjl = V" 2 ? m 
agreement with the conclusion of [Tol |. 



It is easy to see that the reason this above strategy 
does worse than the shot- noise limit (SNL), which can be 
achieved with simple photon intensity difference count- 
ing 0, is that the projective measurement on the opera- 
tor Nab, defined above, throws away very many photon 
number amplitudes in the coherent state, keeping only 
the two terms in Eq. (2) with n — N and m = or n = 
and m = N. Such an observation suggests that we intro- 
duce a new operator, which projects onto all of the maxi- 
mally super-resolving terms, that is, vab = Y^n=o Nab, 
which corresponds to the phase-bearing anti-diagonal 
terms in the two-mode density matrix describing the in- 
terferometer [(| . We may now carry out the expectation 
with respect to the dual-mode coherent state of Eq. (2) 
to get, 

{v A b) = 2e- fi+ " COS¥,/2 cos(nsin^/2) 

^ 2e- fi/2 cos(r^/2), (5) 

where we have approximated the expression near ip — 
0, where the function is sharply peaked. We can see 
from this expression that by carefully choosing the return 
power n = N, we recover super resolution. However 
the exponential pre-factor produces a loss in visibility. 
Choosing n — N = 20, to obtain again 10-fold super- 
resolution, we obtain a visibility of only 0.009%. This 
low visibility does not bode well for the sensitivity, which 
we may compute as per Eq. (3) to get, 

A^ 2 = |e fi + e 3 "/ 2 -e" COS¥) [l + cos(nsin^)]| 

(Vfsin^)A, (6) 
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which again displays the tempting Heisenberg-limiting 
pre-factor of 1/n 2 . This expression is singular at the 
phase origin, but has a minimum nearby ip — 7r/2n, which 
for large photon number is approximately at n = N, 
which by again choosing this to be a large integer, sim- 
plifies Eq. (6) to, 
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which only approaches the SNL and HL near n — 1 and 
then rapidly diverges to be much worse than the SNL 
for large n. Hence, from these examples, we see there is 
a high price to pay for NOON-like super-resolution with 
coherent states — so many photon amplitudes are dis- 
carded that we always do far worse than shot noise. This 
analysis then suggests our final protocol — what if wc 
choose a measurement scheme that includes all of the 
phase-carrying off-diagonal terms in the two-mode den- 
sity matrix? Such a scheme is to consider the operator 
constructed from all the MM' projectors [6(, that is, 
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FIG. 2: This plot shows the expectation value ((lab) of Eq. 
(9) plotted as a function of the phase shift ip (solid curve) for a 
return power of n — 100. For reference we plot the normalized 
"classical" two-port difference signal (dashed curve). We see 
that the plot of the (P>ab) is super-resolving and is narrower 
than the classical curve by a factor of Sip = = 1/10, as 

given in Eq. (11). 
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FIG. 3 
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In this plot we depict the sensitivity expression Atpt 
(10), again for the return power of n = 100 (solid 
The horizontal dashed line indicates the shot-noise 
limit of A(/3g NL = 1/n = 1/100. We see that the sensitivity 
of the super-resolving (lab detection scheme hits the SNL at 
ip = 0, as indicated by expanding Eq. (10) in a power series. 



where we note this is evidently not a resolution of the 
identity operator. It is easy to show that this operator 
of Eq. (8) is both Hermitian and idempotent, that is 
fi AB = fiAB and fX 2 AB = Iab, respectively, where Iab is 
the two-mode identity operator. Using these properties, 
with a bit of algebra, we establish, with respect to the 
two-mode coherent state of Eq. (2), the following results, 
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where once again the Heisenberg limit term of 1/n 2 ap- 
pears accompanied by an exponential factor in the sensi- 
tivity estimate. We plot as a solid curve the expectation 
value of Eq. (9), for a return power of n = 100, in Fig. 2, 
along with the standard (classical) photon difference de- 
tection interferogram (dashed curve) 0]. Clearly {(iab) 
has a visibility of 100% now, and is periodic in tp with 
period 27r, and highly peaked at the phase origin where 
ip = 0. This curve is not super-resolving in the usual 
sense of the word, as there are no multiple narrow peaks 
as would be the case in a NOON-state scheme, but it is 
super-resolving in the sense that there is a well defined 
narrow feature that is clearly sub-Rayleigh limited in res- 
olution. Such a feature would be useful, for example, in 
LADAR ranging or laser Doppler velocimetry, where one 
would lock onto the side of such a feature and then mon- 
itor how it changes in time with a feedback loop in the 
interferometer. 

To estimate the width of this central peak we note that 
in the small phase angle limit, Eq. (9) may be approxi- 
mated as, 
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by choosing a return power of n — 100, we are 10-fold 
super resolving in terms of this central feature of the 
expectation value. 

Now we check the sensitivity of this scheme. In Fig. 3 
we plot as a solid curve the variance of Eq. (10), again 
for n = 100, near the phase origin over — tt/8 < (f < ir/8. 
We include as a dashed curve the shot-noise limit. We 
see that the sensitivity is shot-noise limited about the 
phase origin. This may be established analytically by 
first noting that the expression of Eq. (10) for the sensi- 
tivity has a removable singularity at the origin, and then 
by expanding it in a power series around <p> = to get 
^V/il^so = 1/v^i which is precisely the shot-noise limit. 
Hence by counting all the photons in an off-diagonal fash- 
ion, the detection scheme embodied in the operator (lab 
of Eq. (8) produces a new kind of super resolution, and 
performs at the shot-noise limit in sensitivity. Since all 
of this information is extracted at the detector, and only 
coherent states are used at the source and in the interfer- 
ometer, this scheme will be no worse in sensitivity in the 
presence of absorption or loss than an equivalent classi- 
cal LADAR scheme, but it will in addition have super- 
resolving capabilities. 

It remains to be understood how the observable of Eq. 
(8) might be detected in the laboratory. We note that 
at this stage of the analysis the operator (iab is to be 
carried out with respect to the two-mode coherent state 
at the line III in the Fig. 1. That is we have not yet 
applied the second beam splitter. To the right of the 
second beam splitter, at the line IV in Fig. 1, we wish 



to now carry out the well-known parity operator ll| on 
output mode A in the upper arm, 



n.4 = (-i) 



7iA irra^a 
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which is clearly a Gaussian of width 6p> — 1/yn. Hence 



which simply indicates whether an even or odd number 
of photons exits that port. Here ua = a' a is the num- 
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ber operator for that mode. Such a detection scheme is 
easily implemented by placing a highly efficient photon- 
numbcr-resolving detector at this port, and such detec- 
tors with 95% efficiency and number resolving capabili- 
ties in the tens of photons have been demonstrated 
The connection between these two operators may be es- 
tablished via the easily proved identity, 



AB \ 



I MAS I 



I AB 



(a,0\U\n A ®i B )U\a,0), (13) 



where Ib is the _B-mode identity operator, and U = 
e~ tJ yf with J y — i{ab* — a^b)/2 denotes the transfor- 
mation of the usual Mach-Zehnder interferometer. That 
is, the effect of measuring the two-mode coherent state 
with respect to the operator fiAB to the left of the second 
BS at line III is equivalent to measuring it with respect 
to the parity operator 11^ to the right of the second BS 
at line IV, all indicated in Fig. 1. Hence Eq. (14) estab- 
lishes that the two schemes have the same super-resolving 
and shot-noise limiting properties, with the important 
point that the parity operator is perhaps far easier to 
implement in the laboratory, and it has been show to be 
a universal detection scheme in quantum interferometry 

II, El- 

In summary we have provided a super-resolving inter- 
ferometric metrology strategy, which achieves the shot- 
noise limit. The protocol has the appealing feature that 
it requires only the production and transmission of ordi- 
nary laser beams in the form of coherent states of light. 
Hence, unlike the issues concerning the propagation of 
non-classical states of light, such as squeezed light or en- 
tangled Fock states, this scheme clearly suffers no worse 
degradation in the presence of absorption and loss than 
a classical coherent LADAR system. All of the quantum 
trickery, which provides the super resolution, is carried 
out in the detection, which can be carried out with cur- 
rent photon number-resolving technology. However, a 
scheme by which we count the number of photons and 
then decide if that number is even or odd is overkill. 
The parity operation only requires that we know the sign 
— even or odd — independently of the actual number. 
Hence counting photon number is sufficient but perhaps 
not necessary, if a general scheme to determine the par- 
ity of a photon state could be found that did not require 
photon number counting. We conjecture such a scheme 
exists, perhaps through the exploitation of optical non- 
linearities [15] ], or projective measurements, and this is 
an area of ongoing research. Our protocol for super- 
resolving phase measurements at the shot-noise limit has 
applications to quantum imaging, metrology, and remote 
sensing. In particular, for applications such as Dopplcr 
velocimetry of rapidly moving objects, one does not have 
the luxury to integrate the data for long periods of time 
in order to push down the signal to noise. In such sce- 



narios the signal resolution in the form of the Rayleigh 
criteria is the usual limit to the system performance. It 
is in such situation that we anticipate this protocol to be 
most useful. 
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